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Abstract 


The lacertid genus Pedioplanis is a moderately speciose group of small-bodied, cryptically-colored lizards found in arid habitats 
throughout southern Africa. Previous phylogenetic work on Pedioplanis has determined its placement within the broader context of 
the Lacertidae, but interspecific relations within the genus remain unsettled, particularly within the P. undata species complex, a group 
largely endemic to Namibia. We greatly expanded taxon sampling for members of the P. undata complex and other Pedioplanis, and 
generated molecular sequence data from 1,937 bp of mtDNA (ND2 and cyt 5) and 2,015 bp of nDNA (KIF24, PRLR, RAG-1) which 
were combined with sequences from GenBank resulting in a final dataset of 455 individuals. Both maximum likelihood and Bayesian 
analyses recover similar phylogenetic results and reveal the polyphyly of P. undata and P. inornata as presently construed. We con- 
firm that P. husabensis is sister to the group comprising the P. undata complex plus the Angolan sister species P. huntleyi + P. haackei 
and demonstrate that P. benguelensis lies outside of this clade in its entirety. The complex itself comprises six species including P. 
undata, P. inornata, P. rubens, P. gaerdesi and two previously undescribed entities. Based on divergence date estimates, the P. undata 
species complex began diversifying in the late Miocene (5.3 + 1.6 MYA) with the most recent cladogenetic events dating to the Plio- 
cene (2.6 + 1.0 MYA), making this assemblage relatively young compared to the genus Pedioplanis as a whole, the origin of which 
dates back to the mid-Miocene (13.5 + 1.8 MYA). Using an integrative approach, we here describe Pedioplanis branchi sp. nov. and 
Pedioplanis mayeri sp. nov. representing northern populations previously assigned to P. inornata and P. undata, respectively. These 
entities were first flagged as possible new species by Berger-Dell’mour and Mayer over thirty years ago but were never formally 
described. The new species are supported chiefly by differences in coloration and by unique amino acid substitutions. We provide 
comprehensive maps depicting historical records based on museum specimens plus new records from this study for all members 
of the P. undata complex and P. husabensis. We suggest that climatic oscillations of the Upper Miocene and Pliocene-Pleistocene 
era in concert with the formation of biogeographic barriers have led to population isolation, gene flow restrictions and ultimately 
cladogenesis in the P. undata complex. 
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Introduction 


The Lacertidae is a large family of Old World lizards with a 
center of diversity in southern Africa, particularly in xero- 
thermic habitats (Arnold et al. 2007; Mayer and Pavlicev 
2007; Garcia-Porta et al. 2019). Overall, this family has a 
conservative general morphology and the degree of homo- 
plasy and convergence in external features 1s very high, es- 
pecially in predominantly arid-zone lineages (Borsuk-Bi- 
alynicka et al. 1999; Hipsley and Miller 2017). With 13 
currently recognized species, Pedioplanis Fitzinger, 1843 
is the most diverse genus of lacertid lizards in southern Af- 
rica (Conradie et al. 2012), with taxa occupying semi-arid 
to highly xeric habitats, including coastal fynbos, Nama 
and succulent Karoo, grassland, savannah, and_hard- 
packed gravel plains and deserts (Branch 1998; Makokha 
et al. 2007). Nevertheless, previous studies have detected 
additional, novel lineages among Pedioplanis species that 
await resolution and formal description (Mayer and Berg- 
er-Dell’mour 1987; Makokha et al. 2007; Conradie et al. 
2012). Members of the Pedioplanis undata species com- 
plex are autochthonous to Namibia and South Africa and 
have proven particularly taxonomically challenging (May- 
er and Berger-Dell’mour 1987; Makokha et al. 2007). All 
are similar in size (adult SVL 40-56 mm) and overall 
body form, but exhibit differences in dorsal patterning 
and coloration, arrangement of scales in the lower eyelid, 
and geographic distribution (Mayer and Berger-Dell’mour 
1987; Branch 1998; Makokha et al. 2007). 

Pedioplanis undata was originally described as a single 
widespread species, Lacerta undata Smith, 1838, and one 
year later transferred to the genus Eremias Fitzinger, 1834 
by Dumeéril and Bibron (1839). Three non-nominotypical 
subspecies, Eremias undatarubens, E. u. gaerdesi and E. u. 
inornata (the last originally described as full species, Ere- 
mias inornata Roux, 1907) were subsequently recognized 
based on morphological differences (Boulenger 1921; 
Mertens 1954, 1955). Mayer and Berger-Dell’ mour (1987) 
used allozyme data to delineate a total of seven geographi- 
cally coherent forms including northern and southern vari- 
eties of P. u. undata and P. u. inornata. They provisionally 
elevated P. gaerdesi and P. rubens to full species, and later 
used the same methods to describe P. husabensis, which 
together with P. u. undata and P. u. inornata resulted in 
five putative members of the complex (Berger-Dell’mour 
and Mayer 1989). Morphologically-based phylogenetic 
analyses by Arnold (1991) recovered the same five named 
taxa at the species level, although their results suggested 
that the P. undata complex was rendered paraphyletic by 
P. namaquensis and P. benguelensis. 

Multiple recent molecular phylogenies have con- 
firmed the monophyly of the P. undata species complex 
(including P. husabensis, as sister to all remaining spe- 
cies; Makokha et al. 2007; Conradie et al. 2012), how- 
ever, some interspecific relationships remain unresolved. 
Maximum likelihood analyses place P. rubens in a clade 
with P. undata that is sister to all remaining species with- 
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in the complex (Makokha et al. 2007), whereas Bayesian 
analyses suggest that the former is more deeply nested 
within the complex (Makokha et al. 2007; Conradie et 
al. 2012). The monophyly of P. inornata has also been an 
ongoing subject of debate. Makokha et al. (2007) identi- 
fied two moderately divergent (ND2 = 7-8%; 16S rRNA 
= 3-4%; RAG-1 = 1—2%) lineages within the species, a 
west-central Namibian clade and a southern Namibian/ 
Northern Cape (South Africa) clade. Pedioplanis gaerde- 
si was recovered as the sister to the former, rendering P. 
inornata paraphyletic (Makokha et al. 2007). These re- 
sults corroborate Mayer and Berger-Dell’mour’s (1987) 
report of two morphologically and geographically dis- 
tinct forms of P. inornata, however, no formal taxonomic 
changes were made due to limited geographic sampling. 
Finally, the status of P. undata has also previously come 
into question, as Mayer and Berger-Dell’mour (1987) 
identified northern and southern forms of the species, yet 
the most recent phylogenetic analyses detected no such 
divergent clades given their limited sampling (Makokha 
et al. 2007; Conradie et al. 2012). 

We present the most comprehensive phylogenetic as- 
sessment yet of the P. undata species complex in order to 
resolve interspecific relationships among described species, 
investigate the status of previously recognized unnamed 
lineages, and elucidate the geography of species diversifi- 
cation within the complex. We compiled a new and extend- 
ed multi-locus dataset and reconstructed the phylogeny of 
the genus Pedioplanis as a whole with a greatly increased 
sample set for each member of the P. undata species com- 
plex throughout their ranges, particularly in Namibia. 


Materials and methods 
Field sampling 


The majority of specimens included in this study were 
collected by the authors, under permitted fieldwork in 
Namibia, South Africa and Botswana. Specimens were 
euthanized, fixed in 10% neutral buffered formalin, and 
stored in 70-75% ethanol (EtOH). Prior to fixation, tail, 
thigh muscle or liver tissue was preserved in 95% EtOH 
for subsequent DNA sequencing. Specimens have been 
deposited in the Museum of Comparative Zoology at 
Harvard University (MCZ), the California Academy of 
Sciences (CAS), the National Museum of Namibia in 
Windhoek (NMNW), and the Museum fiir Naturkunde in 
Berlin (MfN, ZMB). All locality data were recorded by 
means of handheld GPS devices in decimal degrees (DD) 
using the WGS 84 reference coordinate system. 


Molecular phylogeny 


Genomic DNA was extracted using the Qiagen DNAeasy 
Kit. PCR amplification was performed on an Eppendorf 
Mastercycler gradient thermocycler using primer pairs 
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obtained from published sources (see Suppl. material 4: 
Table S1 for primer pairs). Molecular sequence data were 
generated for two mitochondrial genes, NADH dehydro- 
genase subunit 2 (ND2) and cytochrome b (cyt b), and 
for three nuclear genes, recombination activation gene 1 
(RAG-1), prolactin receptor (PRLR), and kinesin family 
member 24 (KIF24). Nuclear genes were sequenced for 
only a subsample of individuals representing unique mi- 
tochondrial haplotypes (KIF24: n= 229; PRLR: n= 214; 
RAG-1: n = 244). PCR products were visualized using 
1.5% agarose gels and then purified with a homemade 
magnetic bead solution (Rohland and Reich 2012). Cycle 
Sequencing was performed using the BigDye Terminator 
v3.1 Cycle Sequencing Kit (Applied Biosystems, Foster 
City, CA, USA) and purified with an additional magnetic 
bead cleanup before being sequenced on an ABI 3730xl 
DNA analyzer. Sequences were assembled and edited in 
Geneious v9.1.8; initial sequence alignments were made 
using MUSCLE v3.8.31 (Edgar 2004) before adjusting 
by eye. Sequence alignments were used to calculate av- 
erage uncorrected pairwise genetic distances (p-distance) 
using the distance matrices generated in Geneious v9.1.8. 

Sequences were generated for 428 individuals, includ- 
ing all 13 currently recognized species of Pedioplanis 
and eight outgroup taxa (Nucras aurantiacus, N. holubi, 
N. tessellata, Heliobolus lugubris, Ichnotropis capensis, 
Meroles knoxii, M. reticulatus, M. suborbitalis). Based 
on recent molecular phylogenies (Mayer and Pavlicev 
2007; Engleder et al. 2013) the sister pair Nucras + 
Heliobolus were used to root the tree. Within the P. un- 
data species complex, the number of individuals for 
which sequences were generated are as follows: 18 
P. husabensis (14 localities, see Kirchhof et al. 2017), 
71 P. undata (28 localities), three P. rubens (2 localities), 
110 P inornata (58 localities), 48 P. gaerdesi (37 local- 
ities). Additional Pedioplanis sequences from the genes 
ND2 (n = 45) and RAG-1 (n = 40) were downloaded 
from GenBank, providing thorough representation of 
each species across its known geographic range (see 
Suppl. material 1 for complete list of tissue specimen 
and GenBank sampling). 

Phylogenetic analyses using Maximum likelihood 
(ML) and Bayesian inference (BI) were performed on the 
CIPRES Science Gateway v3.3 (Miller et al. 2010). Sep- 
arate analyses were performed on individual genes and 
on the concatenated sets of mitochondrial and nuclear 
genes separately. Best-fitting models of evolution were 
identified using PartitionFinder v.2.1.1 (Lanfear et al. 
2016) in which partitioning schemes were selected using 
the Bayesian Information Criterion (BIC) score for each 
analysis. All ML partitions were run under the GTR+IT 
model of evolution using RAxML v8.1.24 (Stamatakis 
2014) for 1,000 rapid nonparametric bootstrap replicates; 
BI partitions were run under the models listed in Sup- 
pl. material 5: Table S2. MrBayes v3.2.6 (Ronquist and 
Huelsenbeck 2003) was used to perform all Bayesian 
phylogenetic analyses and run for 50,000,000 genera- 
tions, sampling every 10,000 generations; convergence 
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of the Markov chains was assessed by eye using Tracer 
v1.7 (Rambaut et al. 2018) and the initial 25% of trees 
were discarded as burn-in. 


Morphological data and analysis 


Morphological comparisons among members of the P. un- 
data complex were made following initial phylogenetic 
analysis. A diversity of mensural and meristic measures 
was recorded for 120 specimens including the type series 
of the new species described herein, and additional speci- 
mens representative of major clades revealed by the phy- 
logenetic analyses. In order to obtain adequate sampling 
to capture morphological variation in each clade, we add- 
ed non-genotyped adult individuals (> 40 mm snout—vent 
length) of P. undata sensu stricto (n= 6) collected from the 
Oanob Dam (near Rehoboth) population, and of P. rubens 
(n = 9) from the Klein Waterberg population, to the mor- 
phological analyses (see Appendix | for a complete list of 
specimens used for the morphological analyses). Morpho- 
logical data included the following measurements in mm 
(for bilateral features only right side values were record- 
ed, unless damaged): Snout—vent length (SVL) — snout tip 
to the anterior edge of the cloaca; body length (BL) — ante- 
rior edge of cloaca to posterior edge of collar; collar-snout 
length (CSL) — posterior edge of collar to snout tip; head 
length (HL) — posterior edge of the occipital to snout tip; 
maximum head width (HW); lower jaw length (JL) — an- 
terior edge of jaw bone to tip of lower jaw; length of the 
fourth finger (FiL) and fourth toe (ToL); inter-limb length 
(IL} length between axillary and inguinal regions; hind 
limb (tibia) length (HiL); forearm length (FoL) — elbow 
to wrist; tail length (TaL) — posterior edge of cloaca to tail 
tip. All measurements were collected using Mitutoyo dig- 
ital calipers (to 0.1 mm). Comparative data for Angolan 
taxa were obtained from Conradie et al. (2012). 

Features of head and body scalation recorded were: 
prefrontal contact (PF); number of small granules in 
front of supraoculars touching frontal and prefrontal 
(G); number of rows of granules between supraocu- 
lars and supraciliaries (RG); number of supraciliaries 
(Su); interparietal-occipital contact (IO); number of 
infralabials (IF); number of gular scales in a straight 
line between the chin symphysis and median collar 
plate (Gu); number of subdigital lamellae on fourth toe 
(SuL); number of femoral pores on right leg (left if 
right leg was damaged or pores obscured) (Fe). Fea- 
tures of dorsal color and pattern were also recorded. 
All recorded morphological data for the type spec- 
imens are reported in the species descriptions below 
and listed in Tables 1 and 2; summaries of the meristic 
and mensural character data for members of the P. un- 
data species complex can be found in Suppl. materi- 
als 2, 3 (Suppl. material 2: mensural character data for 
P. undata species complex voucher specimens; Suppl. 
material 3: meristic character data for P. undata spe- 
cies complex voucher specimens). 
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X-rays and Micro-computed (micro-CT) scans of spec- 
imens in the type series designated in this study were pre- 
pared at the California Academy of Sciences and Museum 
ftir Naturkunde, respectively, to obtain presacral vertebrae 
counts (PV). All x-ray images were taken using a Faxitron 
Cabinet X-ray System Model 43855C (Faxitron, Tucson, 
AZ, USA) at 45 kV for 60 s. Micro-CT scans were ob- 
tained using a Phoenix nanotom X-ray|s tube at 85 kV and 
120 or 150 yA, generating 1440 projections with 750 ms 
per scan, with differing kV-settings depending on the 
respective specimen size. Effective voxel size, 1.e., res- 
olution in three-dimensional space, ranged from 12.8 to 
14.3 um. The cone beam reconstruction was performed 
using the datos|x-reconstruction software (GE Sensing & 
Inspection Technologies GMBH phoenix|x-ray datos|x 
2.2) and the data were visualized in VG Studio Max 3.1. 
Micro-CT scans have been deposited on the MfN data re- 
pository (https://doi.naturkundemuseum. berlin/10.7479/ 
4)¢6-ae46) and are available on the MorphoSource web- 
site (http://morphosource.org). X-ray images are also 
available on MorphoSource (See Appendix | for DOIs). 

Due to limited sample sizes we ran our morphological 
analyses on a dataset combining data for adult males and 
females of each species; all statistical analyses were per- 
formed in R (R Core Team 2019). Meristic variables were 
tested for normal distribution using the Shapiro-Wilk nor- 
mality test as part of the R package dplyr (Wickham et al. 
2021) and for homogeneity of variance with Levene’s test 
using the package car (Fox and Weisberg 2019). They were 
then compared between species using one-way ANOVA (for 
normally distributed data with homogeneity of variance) or 
Kruskal-Wallis rank sum test (for skewed data or data with- 
out homogeneity of variance). We excluded the categorical 
variables prefrontal contact (PF) and interparietal-occipital 
contact (IO) and number of rows of granules between su- 
praoculars and supraciliaries (RG) from this analysis. 

Allometric effects were taken into account when men- 
sural characters were compared. We built bivariate linear 
models using the raw measurements of each character in 
relation to the respective individual’s SVL and compared 
the slopes of the regression lines visually using the ggplot2 
package (Wickham 2016). We refrained from additional 
multivariate analyses due to a concern that doing so would 
potentially lead to unreliable results as a consequence of 
comparatively smaller sample sizes in the datasets. 


Results 


Phylogenetic analysis 


Final alignments of the in-group Pedioplanis taxa were 
as follows: ND2: 905 base pairs (429 variable, 379 parsi- 
mony informative); cyt b: 1,032 base pairs (444 variable, 
363 parsimony informative); KIF24: 484 base pairs (72 
variable, 50 parsimony informative); PRLR: 544 base 
pairs (73 variable, 46 parsimony informative); RAG-1: 
987 base pairs (124 variable, 74 parsimony informative). 
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There were no conflicts in the topologies between the 
BI and ML analyses and both retrieved generally high 
nodal support (bootstraps — BS > 70%; posterior prob- 
abilities — PP > 0.95) throughout their respective trees. 
Separate analyses among the nuclear genes revealed no 
conflicts, though individually they provided little resolu- 
tion and were thus concatenated into a single alignment. 
Within the concatenated nuclear tree (not shown) P. ben- 
guelensis, P. husabensis and P. rubens were found to be 
monophyletic with strong support (P. benguelensis: PP = 
1.0, BS = 100%; P. husabensis: PP = 1.0, BS = 100%; 
P. rubens: PP = 0.85, BS = 100%). However, all supra- 
specific clades within the P. undata complex were weakly 
supported by our nuclear data, and relationships among 
them remained largely unresolved. 

Separate analyses were performed on the individual 
and concatenated mitochondrial genes ND2 and cyt 5. 
Each analysis resulted in a well-resolved tree and no tax- 
onomically relevant topological conflicts were revealed. 
Our concatenated mitochondrial analyses (Fig. 1) recov- 
ered a strongly supported monophyletic Pedioplanis as 
sister to Meroles + Ichnotropis (PP = 1.0, BS = 98%). 
Within Pedioplanis, P. burchelli + P. laticeps are sister 
to all remaining congeners (PP = 1.0, BS = 100%). Pe- 
dioplanis breviceps and P. lineoocellata are recovered as 
sister taxa (PP = 1.0, BS = 100%) and are, in turn, sister to 
P. namaquensis plus all remaining Pedioplanis (PP = 1.0, 
BS = 100%). Pedioplanis benguelensis, P. husabensis and 
(P. huntleyi + P. haackei) are successively more closely re- 
lated to the P. undata complex (PP = 1.0, BS = 99-100%). 

Within the clade containing the four recognized species 
of the P. undata species complex (P. undata, P. gaerdesi, 
P. inornata, P. rubens) we recovered a total of six highly 
divergent lineages (Figs 1, 2). In all analyses P. undata 
appears polyphyletic, comprising two geographically 
distinct lineages, one from northern Namibia (P. undata 
“North” sensu Mayer and Berger-Dell’mour 1987) and 
another that occurs in central Namibia from the Khomas 
Highlands, south into parts of the northern Hardap Re- 
gion and west into the Erongo Region (P. undata “South” 
sensu Mayer and Berger-Dell’ mour 1987). The two forms 
were not recovered as sister taxa and are highly divergent 
from one another (average uncorrected p-distance + SD: 
cyt b: 12% + 0.9; ND2: 14% + 4.1). The southern lineage 
is deeply nested within the complex, whereas P. undata 
“North” is sister to all other P. undata complex lineages 
(PP = 1.0, BS = 100%). Pedioplanis rubens was recov- 
ered nested within the P. undata species complex and sis- 
ter to the remaining clades (PP = 1.0, BS = 99%). Individ- 
uals that superficially resemble P. inornata or P. gaerdesi 
from central-west Namibia in the Erongo region form a 
clade, (P. inornata “Central” sensu Makokha et al. 2007) 
and are sister to P. gaerdesi sensu stricto (PP = 1.0, BS = 
100%), collected largely from north of the Ugab River. A 
second P. inornata clade (P. inornata “South”) was repre- 
sented by specimens from as far north as the Kuiseb River 
and south into the Northern Cape Province of South Af- 
rica. These P. inornata lineages are moderately divergent 
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Figure 1. Sampling map and phylogeny of Pedioplanis. (a.) Satellite map of Namibia depicting sampling localities for genotyped 
specimens of P. husabensis and the six major clades within the P. undata species complex recovered by maximum likelihood and 
Bayesian analysis of the concatenated mitochondrial dataset. (b, c.) Maximum likelihood phylogram depicting (b.) the interspecific 
phylogenetic relationships among all Pedioplanis species and (c.) all major clades within the P. undata species complex; colored 
triangles indicate multiple sequenced individuals; typical dorsal patterning for each member of the complex is shown based on pho- 
tographs of live individuals and museum vouchers. Bootstrap values (above) and posterior probabilities (below) are provided for 
each major node. Satellite map of Namibia and surrounding countries created by using the Google Maps terrain map layer, which 


South Africa 


was accessed through the Free and Open Source QGIS (QGIS.org 2021) and its OpenLayers Plugin. 


from each other (average uncorrected p-distance + SD): 
cyt b: 7.7% + 0.5; ND2: 8.8% + 4.3). A sister relationship 
between P. undata “South” and P. inornata “South” is 
recovered in all analyses (PP = 0.99, BS = 91%), although 
without support in the case of nuclear data. A sampling 
map for sequenced tissues of P. husabensis and all clades 
recovered within the P. undata species complex is depict- 
ed in Fig. 1, and within these clades, major geographic 
structuring is depicted in Fig. 2a, b. 


Morphological analysis 


The final morphological dataset (adults only) contained 
13 P. undata “North” (7 males, 6 females), 20 P. inornata 
“Central” (15 males, 5 females), 22 P. inornata “South” 
(15 males, 7 females), 10 P. rubens (4 males, 6 females), 
21 P. gaerdesi (8 males, 13 females), and 11 P unda- 
ta “South” (2 males, 9 females). We ran Shapiro-Wilk 
normality tests on all meristic characters and found that 
only Gu data was normally distributed (W = 0.9838, p = 
0.2782) and met the assumption of homogeneity of vari- 
ance (Levene’s test, F = 1.8392, p = 0.113), whereas all 
other characters were not normally distributed: G (W = 
0.95389, p =0.0018), Su (W = 0.56256, p< 0.001), IF (W 
= 0.74376, p< 0.001), SuL (W = 0.95704, p = 0.0044), Fe 
(W = 0.93945, p < 0.001) (see Suppl. material 6: Table S3 
for detailed results). 

Using ANOVA or Kruskal-Wallis rank sum tests (and 
subsequently applied Wilcoxon rank sum test) we found 


significant differences only in the number of subdigital 
lamellae on the fourth toe (SuL; Kruskal-Wallis chi- 
squared = 11.603, p = 0.0407) and the number of fem- 
oral pores on the leg (Fe; Kruskal-Wallis chi-squared = 
24.992, p < 0.001) between the species of the P. undata 
complex. Since none of these differences contribute 
meaningfully to the diagnosis of the new taxa we pres- 
ent the data in the Suppl. material 7: Table S4 for ref- 
erence, but refrain from overinterpreting them in this 
present study. 

For the difficult to distinguish species pairs P. undata 
“North” and P. undata “South”, as well as P. inornata 
“Central” and P. inornata “South”, we assessed mensu- 
ral characters in an allometric context by comparing their 
linear regressions (see Suppl. material 8: Fig. S1, Suppl. 
material 9: Fig. S2 for linear regression comparisons of 
mensural characters). However, in our dataset, the aver- 
age SVL of adult P. inornata (“Central”) and P. inornata 
(“South”) were significantly different, with P inornata 
(“Central”) being significantly smaller (Wilcoxon rank 
sum test, W = 122, p = 0.01401), preventing a meaning- 
ful comparison of all other size-related body measure- 
ments by this method. While the similar-sized P. undata 
(“North”) and P. undata (“South”) (Wilcoxon rank sum 
test, W = 78, p = 0.728), might differ with regard to fore- 
limb length (FoL) and lower jaw length (JL) (see Suppl. 
material 8: Fig. S1), this may be due to other factors (e.g. 
sampling bias, potential sexual size dimorphism), but lim- 
ited sample sizes prevented the use of multivariate analy- 
ses for clarification. 
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Figure 2. Maximum likelihood phylogeny of the Pedioplanis undata species complex inferred from the concatenation of the two 
mitochondrial genes ND2 and cyt 6. Major geographically-structured clades are depicted for each species: (a.) P. branchi sp. nov. 
(yellow), P. gaerdesi (green), P. mayeri sp. nov. (teal), P. rubens (brown), (b.) P. inornata (blue) and P. undata (orange). Triangles 
indicate multiple sampled specimens. Locality names in bold refer to more inclusive area represented by the clade, parenthetical 
place names are main specific localities. Number ranges correspond to the sample numbers in Suppl. material 1: complete list of of 
tissue specimen and GenBank sampling. Node support labels are represented by bootstrap values (BS), with black nodes denoting 
BS > 90, gray nodes denoting 90 > BS > 70, and white nodes denoting BS < 70; bootstrap values of 100 are not shown. 


Systematics 


Our phylogenetic results support the existence of two 
highly divergent, unnamed entities within the P. undata 
Species complex, each previously recognized as a puta- 
tively distinct, but unnamed lineage of a recognized taxon 
(P. undata and P. inornata). The results of the phylogenetic 
analyses further indicate that these lineages exhibit spe- 
cies-level genetic differentiation and, on this basis, we ele- 
vate the P. undata “North” clade and the P. inornata “Cen- 
tral” clade to full species, which we describe here using 
the general lineage concept of species (de Queiroz 1998). 


Pedioplanis undata “South’/P. undata “North” 


Mayer and Bohme (2000) examined the type materi- 
al of Lacerta undata Smith, 1838 and noted that the 
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rediscovered syntype specimens in the National Muse- 
ums of Scotland collection were, in fact, representatives 
of the taxon long associated with the name Pedioplanis Ii- 
neoocellata pulchella (Gray, 1845) and now, based on the 
findings of Edwards (2013), considered as P. lineoocella- 
ta (Dumeril & Bibron, 1839). They proposed to conserve 
the current usage of names within Pedioplanis by calling 
on the ICZN in Case 3085 to set aside the original type 
designation and recognize the designation of a neotype 
for L. undata. Opinion 1992 of the ICZN (Anonymous 
2002) accepted this solution and thus the neotype was 
fixed as Naturhistorisches Museum Wien (NMW) 31886 
from “near Windhoek, Namibia.” This links the name 
P. undata to the P. undata “South” clade. The assignment 
of the name Pedioplanis undata to the more restricted, 
southern clade is notable in that the predominant con- 
cept of the species based on previous studies of P. undata 
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Figure 222. Continued. 


is largely affiliated with populations and representative 
specimens from the more broadly distributed northern 
clade (Branch 1998; Makokha et al. 2007; Conradie et al. 
2012). In these cases the presence of bold dorsal stripes 
was considered diagnostic for the species, however based 
on our results this character is now largely applicable to 
the unnamed northern clade, and is exhibited inconsis- 
tently and with variation among true P. undata. No other 
nomina are currently within the synonymy of Pedioplanis 
undata and none are available for application to the 
P. undata “North” clade, which we describe as: 


Pedioplanis mayeri sp. nov. 
http://zoobank.org/I EDAC87A-1018-4A70-82BF-B9A 0040D0E86 
Figs 3-5 


Holotype (Fig. 3). * Adult 9; MCZ-R193179 (field num- 
ber MCZ-A28704); Namibia, Erongo Region, Omaru- 
ru District, at Farm Omandumba, ca. 1 km W of house 
(-21.49786, 15.62630, 1238 m above sea level (a.s.1.)); 
collected 15 June 2014 by Jackie L. Childers, Aaron M. 
Bauer, William R. Branch, Matthew Heinicke and Johan 
Marais. Holotype and part of type series to be transferred 
to the National Museum of Namibia. 

Paratypes. n = 8 (all specimens are adults unless oth- 
erwise noted); three 9: CAS 214643, MCZ-R185872, 
ZMB 89349; five 5: CAS 214645, MCZ-R185870, 
MCZ-R190213, MCZ-R193125, ZMB 89350 collected 
from various localities in the Kunene Region of north- 
ern Namibia * CAS 214643; 59 km W of Kamanjab 
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Erongo 
(Gobabeb) (134-135) 


Central 
(Namib Naukluft/Ahomas West/ 
Rehoboth/Okahandja) (136-153) 


P. undata 


Northern Cape - SA 
(Upington/Prieska) (154-156) 


Northern Cape - SA 
(Richtersveld) (157-160) 


Hardap Central A 
(Mariental/Kalkrand) (161-177) 


Central Escarpment 
(Solitaire/Namib Rand Reserve/ 
Helmeringhausen) (178-198) 
Karas Central East 


(Kalahari) (199-202) P inornata 


Northern Cape - Karas (SA/NAM) 
(Richtersveld)/Karas East (203-204) 


Hardap Central B 
(Mariental/Gibeon) (205-207) 


Karas Central South 
(Klein Karas/Fish River Canyon/Aussenkehr) 
(208-215) 


Karas Central 
(Keetmanshoop/Konkiep) (216-220) 


Southern Interior 
(Gibeon/Berseba/Aroab-Kalahari) 
(221-231) 


(-19.64871, 14.33984, 1163 m as.l.); collected 1 June 
2000 by Aaron M. Bauer * CAS 214645; same collec- 
tion data as for proceeding * MCZ-R185870; 35 km S 
of Epupa Falls on Okangwati Rd. (-17.26083, 13.21194, 
1035 m a.s.l.); collected 14 August 2007 by Aaron M. 
Bauer, Johan Marais, Ross A. Sadlier, and Stuart V. Niel- 
sen * MCZ-R185872; N Okangwati on Epupa Falls Rd. 
(-17.40361, 13.15861, 1140 ma.s.l.); collected 14 August 
2007 by Aaron M. Bauer, Johan Marais, Ross A. Sadlier, 
and Stuart V. Nielsen * MCZ-R190213; ca. 4 km from 
Swartbooisdrif towards Epembe (-17.38136, 13.82955, 
836 m a.s.l.); collected 27 November 2011 by Aaron M. 
Bauer * MCZ-R193125; collected at the type locality on 
13 June 2014 * ZMB 89349; subadult 2; along the C35 
road, ca. 2 km N of Fransfontein (-20.19976, 15.01453, 
1097 m as.l.); collected 15 August 2014 by Sebastian 
Kirchhof * ZMB 89350; along the C43 road, 29 km N 
of Palmwag (-19.62736, 13.87391, 1124 m as.l.); col- 
lected 15 August 2014 by Sebastian Kirchhof. Elevation 
data for the holotype and paratypes was obtained using 
the GPS or, if not available, from Google Earth (earth. 
google.com) using georeferenced GPS coordinates from 
the collecting localities. 

Diagnosis. Distinguished from P. /ineoocellata, P. la- 
ticeps and P. burchelli by having 10 longitudinal ventral 
scale rows (vs. 12 or more). It is distinct from P. gaerde- 
si, P. benguelensis, P. husabensis, P. namaquensis and P. 
breviceps in possessing a semi-transparent lower eyelid 
with a brille formed by 2-4 enlarged scales (brille formed 
by a single scale in P. benguelensis and P. gaerdesi, low- 
er eyelid with eight opaque scales in P. husabensis and 
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opaque and scaly in P. breviceps and P. namaquensis). 
Dorsal patterning is characterized by the presence of five, 
usually bold dark brown to black, straight-edged dorsal 
stripes, distinguishing it from P rubens (dorsum and 
tail uniform red-brown to brick red, lacking conspicu- 
ous markings with only a hint of a slightly brighter dor- 
so-lateral line on each side), P. inornata, P. gaerdesi and 
P. branchi sp. nov. (dorsum may be light to dark gray be- 
coming gradually more reddish towards the tail, possess- 
ing dark and/or light beige-yellowish speckling but lack- 
ing distinct longitudinal elements), P. haackei (only three 
dark dorsal stripes), and P. undata (dorsal striping bold 
or not, may be reduced with pale longitudinal elements 
or even a single middorsal stripe restricted to the nape). 
It is further distinguished from P. haackei by typically 
having a smaller number of granules anterior to the first 
supraocular (9-16 vs. 12—32), from P. huntleyi in having 
a larger number of granules anterior to the first supraoc- 
ular (9-16 vs. 7-13), and from P. undata in possessing a 
greater maximum number of granular scales anterior to 
the supraoculars (9-16 in P. mayeri sp. nov. versus 8-13 
in P. undata), and a greater maximum number of femoral 
pores on a single leg (12—16 in P. mayeri sp. nov. versus 
11-14 in P. undata) (see Suppl. material 7: Table S4 and 
Suppl. material 3). Pedioplanis mayeri is similar in most 
characteristics to P. huntleyi, which 1s restricted to Ango- 
la, but the latter species is typified by the restriction of the 
five dark dorsal stripes to the anterior half of the trunk, a 
condition found in only some P. mayeri. No single mor- 
phological character unambiguously distinguishes the al- 
lopatric species P. mayeri, P. undata and P. huntleyi. We 
therefore provide several diagnostic characters based on 
the mitochondrial gene ND2 in order to supplement the 
morphological data. Pedioplanis mayeri may be distin- 
guished from all other Pedioplanis species in being char- 
acterized by having the amino acid histidine instead of 
a tyrosine at base pair 61 due to a codon change at that 
position. It also has the amino acid phenylalanine at base 
pair 223, instead of a threonine (P. benguelensis, P. brev- 
iceps, P. haackei, P. huntleyi, P. husabensis, P. laticeps, P. 
lineoocellata, P. namaquensis), an alanine (P. burchelli), 
or a serine (P. branchi sp. nov., P. gaerdesi, P. inornata, P. 
rubens, P. undata) due to a codon change at that position. 
Finally, it has the amino acid threonine due to a codon 
change at base pair 568, instead of a leucine (P. rubens), 
a valine (P. haackei, P. laticeps, P. undata), or an isole- 
ucine (P. benguelensis, P. branchi sp. nov., P. breviceps, 
P. burchelli, P. gaerdesi, P. husabensis, P. namaquensis), 
P. inornata and P. lineoocellata are polymorphic with 
individuals possessing either valine or isoleucine, and it 
is unknown what the molecular character state is for P. 
huntleyi at this position. 

Description of holotype. Body relatively slender; SVL 
45.5 mm; interlimb distance 20.4 mm; femur 10.0 mm; 
tibia 9.8 mm; humerus 5.5 mm; forearm 5.8 mm; body 
length 28.9 mm from groin to collar; collar-snout length 
16.1 mm; fourth finger length 5.7 mm; fourth toe length 
11.1 mm; head narrow and elongated (head width 52% of 
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head length) with narrow pointed snout (width at rear of 
frontonasal 2.2 mm, width at front of eye 3.8 mm); head 
length 11.0 mm; head width 5.7 mm; lower jaw length 
8.5 mm; eye—ear distance 4.1 mm; eye—nostril distance 
3.8 mm; 1.6 mm between the nostrils; original unregen- 
erated tail 116.7 mm. Rostral semicircular, contacting 
nasals and supranasals; nostril surrounded by nasal, su- 
pranasal, and postnasal; nasals unraised relative to ros- 
tral and frontonasal; postnasal contacts nasal, supranasal, 
frontonasal, anterior loreal, and enters the nostril; nos- 
trils circular; two loreals, anterior half the length of the 
posterior; two preoculars; prefrontals in median contact; 
frontal large, with a narrow posterior projection that con- 
tacts the frontoparietals; frontoparietals in medial con- 
tact; interparietal in contact anteriorly with both fronto- 
parietals, laterally with the parietals, and posteriorly with 
the occipital; occipital partially fused with parietals at 
posterior margin; two supraoculars, both in contact with 
the frontal, preceded anteriorly by a group of 15R/16L 
granules, two granules in contact with prefrontal and two 
in contact with frontal (L); one granule in contact with 
prefrontal and two in contact with frontal (R); two rows 
of granules dividing anterior supraocular from supracil- 
iaries; three small scales between last supraciliary and 
parietal; six supraciliaries on each side, anterior-most 
longest; lower eyelid with transparent brille formed of 
two larger, black-edged scales, with a row of five small- 
er scales below; five supralabials anterior to subocular 
and three supralabials posterior to subocular, on both 
sides; subocular bordering lip, its lower edge narrower 
than its upper; 6R/6L infralabials; first infralabial in con- 
tact with the second infralabial, the first chin shield, and 
the mental; four enlarged pairs of chin shields, the first 
three in median contact and enlarging progressively to 
the posterior; two enlarged temporals; tympanum sunk- 
en; no scales projecting significantly past margin of ear 
opening; enlarged narrow scale at anterodorsal margin 
of tympanum. 28 gular scales in a straight line between 
symphysis of chin shields and median collar plate; col- 
lar free, comprising 11 enlarged plates (median trape- 
zoidal) and extending onto side of neck as a crease that 
terminates midway up the lateral surface; dorsal scales 
small, juxtaposed, granular, without keels, lateral scales 
becoming increasingly larger ventrally; 53 rows of gran- 
ular scales around the midbody; ventral plates in 10 lon- 
gitudinal and 30 transverse rows (from collar to groin); 
ventral scales squarish, subimbricate; single transverse 
row of ventrals across chest just posterior to collar lon- 
ger than broad; 11 enlarged precloacal scales, irregularly 
sized; scales on upper surface of forearm large, smooth, 
and overlapping, without keels; scales on lower surface 
of forearm with a series of enlarged plates, at least twice 
width of scales on upper forearm; scales on upper sur- 
face of tibia rhombic, subimbricate and distinctly keeled, 
below with a series of enlarged plates; scales on upper 
surface of femur keeled with enlarged scales along the 
anterior margin; 15 femoral pores on each side; 27 sub- 
digital lamellae under fourth toe; scales on tail uniform, 
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Figure 3. Pedioplanis mayeri sp. nov. adult female holotype (MCZ-R193179) from Farm Omandumba in the Erongo Region of 
Namibia. (a.) Dorsal and (b.) ventral views of the whole body; scale bars indicate 10 mm; (c.) dorsal, (d.) ventral) and (e.) lateral 


views of the head; scale bar indicates 2.0 mm. 


obliquely rectangular in shape, and strongly keeled; 24 
presacral vertebrae. 

Color in alcohol. Dorsum of crown gray-tan with 
small speckles of varying size on head shields; dorsum 
of snout tinged with reddish-brown. Lateral surface of 
head cream-colored with minute gray speckles on loreal 
scales, supralabials and infralabials. Granules bordering 
lower eyelid window white. Color and pattern of tempo- 
ral region confluent and consistent with those on trunk. 
Trunk dorsum with three bold black stripes each bordered 
on either side by whitish-cream stripes that extend to the 
sacrum, becoming indistinct on the tail; middorsal black 
stripe divided anteriorly by a whitish—cream stripe that 
extends from the occiput to the level of the forelimb in- 
sertion, narrowing from a width of five granular dorsal 
scales to two from anterior to posterior. Below the lat- 
eralmost of these pale stripes, which originates on the 
temporal, lies a thick darkly pigmented line to which is 
fused the coalescence of a series of 9R/8L irregularly 
sized ocelli with pale gray-blue centers comprising 6—14 
granules, each surrounded by a brown to black ring. This 
in turn is bordered ventrally by a diffuse white line which 
begins at the corner of the mouth, passes through the ear 
and continues to the groin. Below this line the flanks 
transition to the immaculate pearly white of the venter. 
Dorsal surfaces of forelimbs mottled brown with margins 
of enlarged scales pale; dorsal surfaces of hind limbs me- 
dium brown, enlarged scales on preaxial surface similar 


to those of forelimb, postaxial surfaces bearing scattered 
cream-colored spots, each 1—3 granular scales in extent; 
palms and soles gray with scattered diffuse pinkish to or- 
ange markings. Tail light brown, densely speckled with 
darker markings which continue the bold lateral black 
stripes of the body dorsum as a pair of narrow dark brown 
stripes running along the medial-most keels of the tail for 
at least three quarters of its length. 

Variation within the type series. Measurements for 
the type series are summarized in Table 1 and mensural 
(n = 13) and meristic (n = 17) character data based on 
measurements of additional museum vouchers can be 
found in Suppl. material 2, 3, respectively. Paratypes 
have similar scalation to the holotype except as follows: 
in MCZ-R185872 the prefrontals do not contact and are 
instead divided by a granule; in MCZ-R193125 the in- 
terparietal and occipital are separated by a granule; in 
all paratype specimens the occipital is not fused with the 
parietals; the number of granules in the group preceding 
the anterior supraocular varies from nine to 15; in CAS 
214645 and ZMB 89349 there is one row of granules di- 
viding the anterior supraocular from the supraciliaries 
(usually two); in CAS 214645 and MCZ-R185870 there 
are 7R and 7L infralabials, respectively (usually six); gu- 
lar scales in a straight line between the symphysis of chin 
shields and the median collar plate varies from 27 to 32; 
midbody scale rows vary from 57 to 64; subdigital lamel- 
lae under the fourth toe range from 25 to 28; unilateral 
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Table 1. Morphological data for Pedioplanis mayeri sp. nov. type series. CAS = California Academy of Sciences, MCZ = Museum 
of Comparative Zoology, ZMB = Museum fiir Naturkunde. Abbreviations for character values are as follows: C = In contact; S = 
Separated by granule; SF = Separated by frontal and frontonasal; INC = Tail incomplete. See Materials and Methods for character 
abbreviations and descriptions; Midbody scale rows (MB) were not recorded for all individuals. 


Type Status Holotype Paratype Paratype Paratype 
Museum No. MCZ-R193179 CAS 214643 MCZ-R 185872 ZMB 89349 
Sex Female Female Female Female 
SVL (mm) 45.5 40.7 47.0 43.3 
TaL (mm) 116.7 NA (INC) 108.5 105.4 
PF C C S C 
G 16 9 11 15 
RG 2 2 2 1 
Su 6 6 6 6 
IO C C C C 
IF 6 6 6 6 
Gu 28 28 29 28 
SuL 27 25 25 26 
Fe 15 13 14 13 
PV 24 25 24 25 
MB - - 61 - 


femoral pores number 13 to 14 (right side only); and pre- 
sacral vertebrae number either 24 or 25. 

Coloration in life (Fig. 4). Dorsal coloration and pat- 
tern varied considerably across examined specimens. 
Dorsum of head gray-tan, sometimes possessing dark 
speckles of varying size; dorsum usually bearing three 
bold black, medium-brown or gray stripes bordered on 
each side by white stripes; dorsal stripes often terminate 
at the tail base or begin to fade into a gray and reddish hue 
midway towards the tail base; dorsal stripes may be ab- 
sent entirely, with individuals possessing a uniformly gray 
dorsum, sometimes with dark speckling. All individuals 
examined possess yellow-centered ocelli or yellow spots 
along the flanks, usually bordered ventrally by a white 
stripe and dorsally by a black stripe (or with the dark rims 
of the ocelli confluent with such a stripe); lateral spots 
or ocelli vary in size from ca. 0.5 to 1.0 mm in diameter, 
in some cases smaller spots may merge with the white 
ventral stripe; anterior surface of the thighs and outermost 
ventral plates sometimes tinged or speckled with orange; 
limbs and feet more or less uniformly grayish (front), red- 
dish-brown to brown, or speckled with light or black dots 
or blotches that may form bars; tail usually bright brick or 
orange-red, in some cases dark speckling may be present 
on the dorsal surface of the tail. 

Coloration in alcohol. Pattern as described above 
though dark stripes may fade to dark gray, yellow lateral 
spots may appear cream white or gray-blue, and red col- 
oration usually present on the dorsum of the midbody and 
tail may become a faded orange-pink or may fade to a 
faint, golden-brown hue. 

Distribution. Pedioplanis mayeri sp. nov. is endemic 
to northern Namibia and occurs from south of the Kunene 
River and east of the Namib Desert along the eastern side 
of the escarpment, thence throughout the eastern Kunene 
Region, entering the northeastern parts of the Erongo Re- 
gion and east through the Otjozondjupa Region, reach- 
ing at least as far east as Oshikango (TM17028) in the 
north, Gobabis (Omaheke Region) in the south-east, and 
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Paratype Paratype Paratype Paratype Paratype 
CAS 214645 MCZ-R 185870 MCZ-R 190213 MCZ-R 193125 ZMB 89350 
Male Male Male Male Male 
44.0 51.4 49.2 43.5 52.7 
107.3 124.9 131.8 96.4 121.5 
C C C C C 
14 10 11 10 13 
1 2 2 2 2 
6 6 6 6 6 
C C C S Cc 

6L/7R 7TL/6R 6 6 6 
31 27 29 31 29 
27 28 27 26 26 
14 14 14 14 13 
24 24 24 25 24 
57 64 64 58 - 


Nauchas in the south. It does not enter the Kalahari dune 
fields, and is possibly absent from the Khomas Hochland, 
where it is replaced by P. undata (see Fig. 6 for map of 
locality records). 

Etymology. The specific epithet is a patronym formed 
in the genitive singular honoring our friend and col- 
league, the Austrian lacertid specialist Werner Mayer 
(1943-2015), who first recognized the distinctiveness 
of his namesake species and whose contributions to the 
study of Pedioplanis have been seminal. 

Natural history. Pedioplanis mayeri is widespread 
in northern Namibia and inhabits Tropical & Subtropi- 
cal Grasslands, Savannas & Shrublands and Deserts & 
Xeric Shrublands biomes (following the classification of 
the WWE; Olson et al. 2001), but it is absent from true 
desert. Within the occupied biomes, P. mayeri inhabits 
Namibian savanna woodlands, Kalahari xeric savan- 
na, Kalahari Acacia-Baikiaea woodlands and Angolan 
Mopane woodlands (terminology follows Olson et al. 
2001). Here, the species is active on sandy clay soils 
and more compact, hard soils, rocky flats and slopes, 
in dense grassy patches as well as in more open areas 
interspersed with shrubs (see Fig. 5 for images of char- 
acteristic habitats). In rocky areas and stony slopes indi- 
viduals were found sleeping under stones. The majority 
of the species’ range comprises Namibia’s arid regions 
with a median annual precipitation of 279 mm (61-559 
mm) and a median annual temperature average of 21.1 
°C (18—23.5 °C) (all climate data taken from Fick and 
Hiymans 2017). Pedioplanis mayeri mainly occurs along 
the northern parts of the Great Escarpment at a median 
elevation of 1160 m a.s.l., but reaches altitudes as high 
as ca. 1600 m a.s.l. in Otjozondjupa, and probably as 
low as ca. 200 m a.s.l. in the north along the Kunene 
River. The largest recorded individual was a female 
with SVL = 57 mm (see also Kirchhof et al. 2014; the 
individuals referred to therein as P. undata are now P. 
mayeri). Breeding season appears to be in spring, sim- 
ilar to other southern African Pedioplanis spp. Gravid 
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Figure 4. Life photographs of representative individuals of Pedioplanis mayeri sp. nov. highlighting color pattern variation within 
the species. (a.) Dorsal and (b.) lateral whole body images of an adult possessing bold dorsal stripes collected at the type locality 
(Farm Omandumba, Erongo Region, Namibia), note the row of yellow spots along the flanks; (c.) Dorsal image of an adult from 
the Kamanyab area depicting fainter, medium-brown dorsal stripes; (d.) An adult female collected from Gobabis (ZMB 80391; see 
also Fig. 5a) with dorsal striping of varying boldness and a more grayish-brown hindbody compared to the more reddish hindbody 


observed 1n other individuals. 
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Figure 5. Habitat of Pedioplanis mayeri sp. nov. (a.) Collection site near Gobabis for ZMB 80389-80392 in open savannah wood- 
land, with tall, dense grasses and introduced Opuntia. (b.) Collection site for ZMB 89350 (29 km N of Palmwag, Kunene Region, 
Namibia) in stony Pro-Namib habitat with sparse vegetation dominated by Euphorbia damarana (center). 


females with up to four eggs were found in December 
and February, and small juveniles (SVL = 27 mm) ap- 
peared in mid-February (Kirchhof et al. 2014). Sym- 
patric ground-living diurnal lizards include Trachylepis 
damarana, Meroles squamulosus, Heliobolus lugu- 
bris, Gerrhosaurus flavigularis, and Agama anchietae, 
among others. 


Conservation. The species occurs across a wide area 
of Namibia in which potential threats from agriculture 
and mining are scattered and localized. Local populations 
occur in protected areas within Etosha National Park and 
Waterberg Plateau Park as well as in many local conser- 
vancies. Applying IUCN criteria, we consider P. mayeri 
to be Least Concern. 
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Figure 6. Topographic maps showing locality records for all the species of the Pedioplanis undata species complex and P. husa- 
bensis: (a.) P. husabensis (b.) P. mayeri sp. nov. (¢.) P. rubens (d.) P. branchi sp. nov. (e.) P. gaerdesi (f.) P. undata (g.) P. inornata. 
Black dots represent museum vouchers from various museum collections, those with a white center were sequenced and included 
in the phylogenetic analyses. Museum voucher information was collected from the California Academy of Sciences San Francisco/ 
USA, Carnegie Museum of Natural History Pittsburgh/USA, Ditsong National Museum Pretoria/South Africa, Iziko South African 
Museum Cape Town/South Africa, Museum fiir Naturkunde Berlin/Germany, Museum of Comparative Zoology Harvard Universi- 
ty/USA, Museum of Vertebrate Zoology Berkeley/USA, Natural History Museum of Los Angeles County/USA, Bayworld Museum 
Port Elizabeth/South Africa, National Museum of Namibia Windhoek/Namibia, Zoological Research Museum Alexander Koenig 
Bonn/Germany, in parts via the Global Biodiversity Information Facility (GBIF; https://www.gbif.org/). 
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Pedioplanis inornata “South’/P. inornata 
“Central” 


The type series of Eremias inornata Roux, 1907 comprises 
eight specimens collected from “Oranje-FluB, Kl.-Nam- 
aqualand” [the Orange River, Little Namaqualand, North- 
ern Cape, Province, South Africa]. Syntype ZMA11049 
(Zoological Museum Amsterdam) was subsequently des- 
ignated as the lectotype (Daan and Hillenius 1966). The 
ZMA collections have subsequently been incorporated 
into those of the Naturalis Biodiversity Center in Leiden 
(RMNH). The geographic origin of the specimens and the 
relatively detailed description and accompanying plate, 
clearly tie the name-bearing lectotype of Pedioplanis in- 
ornata to the P. inornata “South” clade. No other nomina 
are currently within the synonymy of Pedioplanis inorna- 
ta and none are available for application to the P. inornata 
“Central” clade, which is here described as: 


Pedioplanis branchi sp. nov. 
http://zoobank.org/7FDBC 1 FF-BEB4-464A-88C5-6CDD19EDS0F2 
Figs 7-9 


Holotype (Fig. 7). Adult 3; CAS 214788 (field number 
AMB 6551); Namibia, S of Karibib at Junction of Road 
D1914 and Road D1952, Karibib District, Erongo Region 
(-22.27038, 15.57471, 1080 m a.s.l.); collected 8 June 
2000 by Aaron M. Bauer. Holotype and part of type se- 
ries to be transferred to the National Museum of Namibia. 

Paratypes. n = 7 (adults); (two 2: ZMB 89310, ZMB 
89311; five J: CAS 214790, CAS 214792, CAS 214794, 
ZMB 89305, ZMB 89316) collected from various local- 
ities located in the Erongo Region of northwestern Na- 
mibia * CAS 214790; S of Karibib at the junction of Rd. 
D19414 and Rd. D1952 (-22.27281, 15.57815, 1075 m 
a.s.l.); collected 8 June 2000 by Aaron M. Bauer * CAS 
214792; same collection data as for proceeding « CAS 
214794; same collection data as for proceeding * ZMB 
89305; 22 km SW of Uis along C35 (-21.42419, 14.76271, 
831 ma.s.l.); 10 August 2014 * ZMB 89316; 22 km SW 
of Uis along C35 (-21.34513, 14.75676, 884 m a.s.l.); 
23 January 2013 * ZMB 89310; Farm Friedhelm Sack 
(-22.52699, 15.54487, 709 m a.s.l.); 14 October 2014; 
¢ ZMB 89311; same collection data as for proceeding; 
all ZMB specimens collected by Sebastian Kirchhof. EI- 
evation data for the holotype and paratypes was obtained 
using the GPS, or, if not available, from Google Earth 
(earth.google.com) using georeferenced GPS coordinates 
from the collecting localities. 

Diagnosis. Distinguished from P. /ineoocellata, 
P. laticeps and P. burchelli by having 10 longitudinal 
ventral scale rows (vs. 12 or more). It is distinct from 
P. benguelensis, P. gaerdesi, P. breviceps, P. namaquen- 
sis and P. husabensis in usually possessing a semi-trans- 
parent lower eyelid with a brille formed by 2—4 scales 
(brille formed by a single scale in P benguelensis and 
P. gaerdesi, lower eyelid with eight opaque scales in P. 
husabensis and opaque and scaly in P. breviceps and 
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P. namaquensis);, in some rare cases P. branchi sp. nov. 
may possess a single transparent scale in lower eyelid, 
those individuals can be distinguished from P. benguel- 
ensis and P. gaerdesi by color and dorsal patterning (see 
below). Dorsal patterning is characterized as being uni- 
formly gray from the mid-back towards the head with 
a reddish hindbody (posterior half of body) and with a 
series of pale to bright yellow spots or ocelli on lower 
flanks, distinguishing it from P. rubens (dorsum and tail 
uniform red-brown to brick red, lacking conspicuous 
markings with only a hint of a slightly brighter dorso-lat- 
eral line on each side), P. mayeri sp. nov., P. haackei, 
P. huntleyi, P. undata (dorsum contains bold stripes or 
other longitudinal elements), P. gaerdesi (never with lat- 
eral ocelli, and speckled with very small black or light 
dots) and P. inornata (spots on flanks are typically pale 
green, not yellow). It can further be distinguished from 
all other Pedioplanis (except P. inornata) in typically 
having a pair of distinct dark markings on the face, one 
through the eye and extending onto the supralabials di- 
rectly below and one more posterior, near the corner of 
the mouth. The new species is significantly smaller than 
P. inornata (P. branchi mean adult SVL = 44.4 mm, max. 
49.1 mm, versus P. inornata mean adult SVL 47.3 mm, 
max. 54.0 mm for specimens sampled here; to 56.0 mm 
elsewhere [Bauer and Shea 2006]). The maturity of 
specimens was confirmed by their mostly or complete- 
ly fused long bone epiphyses. We also provide several 
diagnostic characters based on the mitochondrial gene 
ND2. Pedioplanis branchi sp. nov. can be distinguished 
from all other members of the P. undata species complex 
except P. laticeps and P. undata in being characterized 
by having the amino acid histidine instead of tyrosine 
at base pair 359 due to a codon change at that position. 
It can also be distinguished from all other Pedioplanis 
species except P. Jaticeps in possessing a thymine at 
base pair 300, rather than an adenine (P. benguelensis, 
P. burchelli, P. gaerdesi, P. haackei, P. huntleyi, P. hus- 
abensis, P. inornata, P. mayeri sp. nov., P. namaquen- 
sis, P. rubens, P. undata) or a cytosine (P. breviceps and 
P. lineoocellata). 

Description of holotype. Body relatively slender (SVL 
47.5 mm); interlimb distance 16.8 mm; femur 9.5 mm; 
tibia 8.7 mm; humerus 5.4 mm; forearm 5.7 mm; body 
length 26.9 mm from groin to collar; collar-snout length 
19.3 mm; fourth finger length 4.7 mm; fourth toe length 
10.4 mm; head narrow and elongated (head width 55% 
of head length) with slight constriction at base of rostrum 
and narrow, pointed snout (width at rear of frontona- 
sal 2.2 mm, width at front of eye 4.4 mm); head length 
11.8 mm; head width 6.6 mm; lower jaw length 10.2 mm; 
eye—ear distance 4.4 mm; eye-nostril distance 4.0 mm; 
1.1 mm between the nostrils; complete original tail 
92.3 mm, with epidermis missing from posterior portions. 
Rostral semicircular; contacting nasals; nasals in contact 
medially, unraised relative to rostral and frontonasal; 
postnasal contacts nasal, frontonasal, anterior loreal, and 
enters the nostril; nostrils circular in shape; two loreals, 
anterior loreal half the length of the posterior loreal; two 
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Figure 7. Pedioplanis branchi sp. nov., adult male holotype (CAS 214788) from south of Karibib in the Erongo Region of Namibia. 
(a.) Dorsal and (b.) ventral views of the whole body; scale bars indicate 10 mm; (c.) dorsal, (d.) ventral) and (e.) lateral views of 
the head; scale bar indicates 2.0 mm. 


preoculars; prefrontals in median contact; frontal large, 
with a narrow posterior projection, bordered anteriorly 
by the prefrontals and posteriorly by the frontoparietals; 
frontoparietals in broad medial contact; interparietal in 
contact anteriorly with both frontoparietals, laterally with 
the parietals, and posteriorly with the occipital; occipi- 
tal trapezoidal in shape; two supraoculars, both in medi- 
al contact with the frontal and frontoparietals, preceded 
anteriorly by 8R/8L granules (on left side two granules 
in contact with prefrontal and one in contact with fron- 
tal; on right side two granules in contact with prefrontal 
and three in contact with frontal); single row of granules 
dividing anterior supraocular from supraciliaries, and 
double row of granules dividing posterior supraocular 
from supraciliaries; two small scales between last supra- 
ciliary and parietal; six supraciliaries on each side, the 
anteriormost longest; lower eyelid with transparent brille 
formed of two larger, black-edged scales, with a row of 
five smaller scales beneath; five supralabials anterior to 
subocular and three supralabials posterior to subocular, 
on both sides; subocular bordering the lip, its lower edge 
shorter than its upper; 6R/6L infralabials; first infralabial 
in contact with the second infralabial, the first chin shield, 
and the mental; four enlarged pairs of chin shields, with 
the first three in medial contact and the posterior-most 
largest; no enlarged temporals; tympanum sunk; no scales 
projecting significantly past margin of ear opening; en- 
larged narrow scale at anterodorsal margin of ear open- 
ing. 30 gular scales in a straight line between symphysis 
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of chin shields and median collar plate; collar free, com- 
prising nine enlarged plates (median kite-shaped, pro- 
jecting posteriorly) and extending onto side of neck as a 
crease that terminates midway up the lateral side; dorsal 
scales small, juxtaposed, granular, without keels, lateral 
scales larger towards ventrals; 49 rows of granular scales 
around the midbody; ventral plates in 10 longitudinal 
and 31 transverse rows (from collar to groin); plates of 
the outermost rows squarish, ventral rows usually twice 
as wide as long; single transverse row of ventrals across 
chest just behind collar longer than broad; nine enlarged 
precloacal scales, irregular in shape, median ones larger; 
centralmost enlarged precloacal in contact with six other 
enlarged precloacals; scales on upper surface of forearm 
large, smooth, and overlapping, without keels; scales on 
lower surface of forearm with a series of enlarged plates, 
at least twice the width of scales on upper forearm; scales 
on upper surface of tibia rhombic, subimbricate and dis- 
tinctly keeled, below with a series of enlarged plates; 
scales on upper surface of femur granular and smooth 
with enlarged scales along the anterior margin; 12R/12L 
femoral pores; subdigital lamellae under fourth toe 
26R/28L; scales on tail obliquely rectangular in shape, 
and strongly keeled. 

Color in alcohol. Dorsum of head gray-tan with dis- 
tinct black speckles of varying size on the head shields; 
lateral surface of head with three bold, black vertical 
bars, one at the anterior margin of eye, one originating on 
the ventral eyelid and one just posterior to the orbit, all 
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Table 2. Morphological data for Pedioplanis branchi sp. nov. type specimens. CAS = California Academy of Sciences, MCZ = 
Museum of Comparative Zoology, ZMB = Museum fiir Naturkunde. Abbreviations for character values are as follows: C = In con- 
tact; S = Separated by granule; SF = Separated by frontal and frontonasal; RGN = Tail regenerated. See Materials and Methods for 


character abbreviations and descriptions. 


Type Status Holotype Paratype Paratype Paratype Paratype Paratype Paratype Paratype 
Museum No. CAS 214788 ZMB89310 ZMB89311 CAS 214790 CAS 214792 CAS 214794 ZMB 89305 ZMB 89316 
Sex Male Female Female Male Male Male Male Male 
SVL (mm) 46.6 48 42.5 50.0 40.0 49.0 44.1 46.6 
TaL (mm) 92.3 88.4 111.2 107.0 99.0 99.0 94.1 60.9 (RGN) 
PF C SF C C C € SF 
G 8 12 14 11 11 14 11 9 
RG 1 2 2 2 2 2 1 
Su 6 5 6 6 6 7 zi 6 
IO C C C C C C C Gc 
IF 7 6 6 7 6 6 7 6 
Gu 30 30 26 30 26 30 32 28 
SuL 26 29 28 26 25 25 27 29 
Fe 11 10 12 12 10 13 13 11 
PV 24 25 24 24 24 24 23 24 


extending onto the supralabial scales; body dorsum uni- 
form gray-tan with small, dense, indistinct black speckles. 
The flanks bear reticulate patterning and possess a single 
row of eight ocelli with light gray-blue spots of varying 
size (comprising 8-11 granules) bordered by indistinct 
black rings of 1-3 granules in width; black rings of ad- 
jacent ocelli either separated by 2-4 granules, or in some 
cases connected; venter cream; limbs pigmented above 
similar to dorsum, and white below; dorsum of tail with 
black speckling that is denser than that of the dorsum of 
body; tail transitions to pinkish-yellow towards the end; 
ventral surface of the tail is uniformly light cream. 

Variation within the type series. Measurements for 
the type series are summarized in Table 2 and mensural 
(n = 20) and meristic (n = 33) character data based on 
measurements of additional museum vouchers can be 
found in Suppl. material 2, 3, respectively. Paratypes 
have similar scalation to holotype except: the prefrontals 
are separated by frontal and frontonasal in ZMB 89310 
and ZMB 89316; the number of granules in group pre- 
ceding the anterior supraocular varies from nine to 14; in 
ZMB 89316 there is only one row of granules dividing 
the supraoculars from the supraciliaries; in CAS 214794 
and ZMB 89305 there are seven supraciliaries and five 
in ZMB 89310 on both sides of the head (usually six); 
there are usually 6 to 7 infralabials anterior to the sub- 
ocular; gular scales in a straight line between symphysis 
of chin shields and median collar plate 26—32; subdigital 
lamellae under fourth toe 25—29; femoral pores 10—13 
(RL only); presacral vertebrae 23-25. 

Coloration in life (Fig. 8). Dorsal coloration and pat- 
tern varied considerably across examined specimens. 
Dorsum of head ranges from uniformly gray to distinctly 
speckled with black, sometimes with bold dark blotches 
or vertical bars. Body dorsum gray anteriorly and red-or- 
ange posteriorly, sometimes bearing light or dark spots 
that can be arranged in longitudinal rows, or are irregu- 
larly scattered, or may form horizontal bars or chevrons 
across the posterior dorsum near the tail base; belly uni- 
form cream. The flanks may be patterned with light or 
dark irregular markings, sometimes forming vertical bars 


or chevrons, and are always lined with a more or less dis- 
tinct longitudinal row of yellow spots or yellow-centered 
ocelli bordered by a broken ring of black granules, which 
range from 0.5—1.0 mm in diameter. The anterior surface 
of the thighs sometimes tinged with yellow; limbs and 
feet uniformly gray (front) or reddish-brown (hind), or 
speckled with light or black dots or blotches that may 
form bars. Dorsum of tail is orange-red, often with dense, 
black speckles. 

Coloration (in alcohol). In comparison to the life 
coloration, the row of spots along the flanks that appear 
yellow in life appear gray-blue in preserved specimens 
and may be conspicuous or faint, and the reddish-brown 
color of the posterior body and tail appears darker and 
more grayish. 

Distribution. Pedioplanis branchi sp. nov. is endemic 
to the Erongo Region in Namibia. Its range stretches from 
just south of the Swakop River in the South, where it oc- 
curs in parapatry with P. husabensis, through the pro-Na- 
mib, to the Ugab River and the Brandberg in the north 
and Mount Erongo and Otjimbingwe in the east, probably 
occurring in parapatry with P. mayeri all along its north- 
eastern border (see Fig. 6 for map of locality records). 

Etymology. The specific epithet is a patronym formed 
in the genitive singular honoring our friend and col- 
league, the British-born South African herpetologist, 
William Roy Branch (1946-2018), in recognition of his 
many contributions to African herpetology and in remem- 
brance of many happy trips in the field together. 

Natural history. Pedioplanis branchi inhabits Na- 
mibian savanna woodlands and the gravel plains and 
rocky outcrops in the Namib Desert (Deserts & Xeric 
Shrublands biome; Olson et al. 2001). The species is ac- 
tive on harder substrates, gravel and bare rock (including 
marble, granite, basalt), mostly with only sparse vegeta- 
tion cover and only rarely forages on softer sandy sub- 
strates (e.g. in small washes on slopes) or in more grassy 
patches (Fig. 9). Median annual precipitation within the 
distribution range of P. branchi reaches a maximum of 
only 187 mm (minimum = 16 mm, median = 142 mm), 
while mean annual temperatures are similar to those in 
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Figure 8. Life photographs of representative individuals of Pedioplanis branchi sp. nov. highlighting color pattern variation within the 
species. (a.) Lateral and (b.) dorsal whole body images of an adult specimen from the Chuos Mountains (Chuosberg) (SK376.2014; 
not included in study), note the dense, dark speckling on the dorsum, and bold dark lateral markings. (¢.) Dorsal and (d.) lateral whole 
body images of an adult specimen (SK402.2014; not included in study) from the northward extension of the Swakop River Canyon, 
note the more uniformly gray and red dorsum lacking speckling and fainter yellow spots along the flanks. Lateral (e.) head and (f.) 
whole body images of an adult specimen from the northward extension of the Swakop River Canyon; note the black, vertical bar be- 
neath the lower eyelid, the bold, dark lateral markings on the body, and the distinct longitudinal row of yellow spots along the flanks. 


Figure 9. Habitat of Pedioplanis branchi sp. nov. (a.) Collection site for ZMB 89310 (Farm Friedhelm Sack, Erongo Region, Na- 
mibia) from the extended Swakop River Canyon on a marble/calc silicate slope with sparse grass cover and isolated shrubs. (b.) 


Collection site for ZMB 89612 (not included in study) from a small hill slope in the Tsiseb conservancy in the Pro-Namib, 20 km 
south of the Brandberg (visible in the background) depicting a sparsely vegetated scree slope with a Commiphora shrub. 


i 


the distribution of P mayeri (median = 21 °C, range =50 mm). Breeding season appears to be similar to other 
19.4—22.5 °C). The species occupies lower elevations, southern African Pedioplanis spp. in spring, gravid fe- 
from 77 maz.s.l. near the coast toa maximum of ca. 1190 males were found in December. Sympatric diurnal liz- 
m a.s.l. near the southeasternmost corner of the distribu- ards include, Trachylepis sulcata, Rhoptropus boultoni, 
tion along the Swakop River Canyon (median = 968 m Agama planiceps, Matobosaurus validus, and Varanus 
a.s.l.). The largest recorded individual was a male (SVL = albigularis, among others. 


zse.pensoft.net 


Zoosyst. Evol. 97 (1) 2021, 249-272 


Conservation status. Although P branchi has an 
extent of occurrence below 20,000 km? it occurs in an 
area with low human population density, generally low 
intensity land use, and a high proportion of suitable 
habitat. Populations occur in protected areas within Na- 
mib-Naukluft Park and Tsiseb Conservancy, and possibly 
Dorob National Park. Applying IUCN criteria we consid- 
er P. branchi to be Least Concern. 


Discussion 


Pedioplanis phylogeny and the recognition of 
two new species within the P. undata species 
complex 


In order to resolve the systematics and biogeography of 
the Pedioplanis undata species complex we analyzed a 
new and greatly extended multi-locus molecular data- 
set comprising 455 samples, from all 13 currently rec- 
ognized species of Pedioplanis and eight outgroup taxa. 
The results from our phylogenetic analyses show that in- 
terspecific relationships within the genus Pedioplanis are 
largely congruent with previous studies, but with several 
notable exceptions. We recover with strong support the 
phylogenetic affinities of three species that have hereto- 
fore remained problematic, P. benguelensis, P. husabensis 
and P. rubens. Within the P. undata species complex our 
results support the recognition of two highly divergent 
lineages and multiple instances of polyphyly. 

Historically, P. benguelensis was thought to be closely 
related to P. namaquensis due to similarities in dorsal pat- 
terning (Bocage 1867; Laurent 1964; Makokha etal. 2007). 
More recently, Conradie et al. (2012) recovered P. ben- 
guelensis as sister to the Angolan species pair P. huntleyi + 
P. haackei and the P. undata complex, with P. husabensis 
falling outside of the entire grouping, however these results 
did not receive high nodal support in their analyses. Our 
results confirm the sister relationship between the Angolan 
taxa and members of the P. undata complex, however we 
recover P. husabensis and P. benguelensis as sequentially 
more distant sister taxa to this clade (Fig. 1b, c). Previous 
phylogenetic studies have obtained conflicting topologies 
for the placement of P. rubens. Maximum likelihood anal- 
yses by Makokha et al. (2007) found that P. rubens forms 
a Clade with P. undata (actually, P mayeri sp. nov.) that is 
sister to P. gaerdesi + P. inornata. However, their BI anal- 
ysis recovered P. rubens as nested within the P. undata 
complex, a result corroborated by Conradie et al. (2012) 
and by this present study (Figs Ic, 2a). 

We identify two instances of polyphyly within the 
P. undata complex. First, P. inornata as previously de- 
fined comprises two divergent lineages, one that occurs 
in the eastern and southern parts of Namibia, and in the 
Northern Cape of South Africa (= P. inornata) and an- 
other comprising individuals from central-western Na- 
mibia (= P. branchi sp. nov.). The two lineages are not 
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each other’s closest relatives; rather, P. branchi sp. nov. 
is sister to P. gaerdesi. This result was corroborated by 
the nuclear data, in which the concatenated phylogeny 
recovered these clades as being distinct from one anoth- 
er. Interestingly, some P. branchi sp. nov. individuals 
possess a single (instead of two) large semi-transparent 
scale (“brille”) in the lower eyelid, which caused them 
to be identified as P. gaerdesi in the field given that this 
scale character 1s currently used to diagnose the species 
(Conradie et al. 2012). We suggest that this morpholog- 
ical character alone is inadequate for diagnosing species 
within the P. undata complex. Makokha et al. (2007) 
identified the presence of a P. inornata “Central” clade 
(= P. branchi sp. nov.), however, due to limited sam- 
pling they were unable to identify its geographic extent. 
Here we provide evidence that this clade 1s limited to an 
arid region north of the Kuiseb River and south of the 
Ugab River. These results corroborate Mayer and Berg- 
er-Dell’mour’s (1987) hypothesis that within P. inornata 
there exist two geographically and morphologically dis- 
tinct forms, one restricted to west-central Namibia that 
possesses a gray forebody and red-orange hindbody with 
yellow lateral spots (P. inornata-“Central” = P. branchi 
sp. nov.; body coloration depicted in Fig. 8), and a more 
broadly distributed form in southern Namibia that is char- 
acterized by brownish or reddish coloration with greenish 
lateral spots, P. inornata-“South” (= P. inornata). 

Another equally striking result is the discovery that 
P. undata comprises two genetically distinct and polyphy- 
letic lineages. The first clade, P. undata “North” (= P. may- 
eri Sp. nov.) includes individuals primarily from northern 
Namibia, (although a single sample was found as far south 
as Nauchas), and is sister to all remaining members of the 
P. undata species complex. The second lineage, P. undata 
“South” (= P. undata sensu stricto), is more deeply nested 
within the complex and is closely related to P. inornata, 
a result recovered by both the nuclear and mtDNA data. 


Diversification and geographic context 


Geographic isolation of populations through vicariance 
can result from various events, including both geologi- 
cal (formation of rivers, valleys, mountains, dunes) and 
climatic (drastic environmental change resulting from 
cooling, warming, or changes in precipitation resulting 
in habitat expansion or contraction). Southern Africa has 
experienced a complex paleoclimatic and geological past. 
In particular, the climatic oscillations of the Upper Mio- 
cene and the Pliocene-Pleistocene era are proposed to 
have been influential in the evolution and diversification 
of southern African fauna (e.g. Bauer 1999; Matthee and 
Flemming 2002; Montgelard and Matthee 2012). 

In addition, multiple potential geographical barriers 
have been hypothesized to be responsible for vicariance 
and subsequent cladogenesis. These include the Great 
Escarpment, the sand seas of the Namib Desert, shifting 
Kalahari sands, and rivers such as the Kuiseb, Zambezi, 
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and Orange (e.g. Lamb and Bauer 2000; Matthee and 
Flemming 2002; Scott et al. 2004; Bauer et al. 2006). 

The current distribution ranges of members of the 
P. undata complex (see Fig. 6) do not indicate that the 
aforementioned barriers today have a substantial effect in 
restricting gene flow. For example, the sister taxa P. inor- 
nata and P. undata both occur east and west of the escarp- 
ment and instead possess a north-south border (P. undata 
also occurs both north and south of the Kuiseb), and the 
sister pair P. branchi and P. gaerdesi both occur west of 
the escarpment with a north-south border not obviously 
coincident with any geographical barrier. Nevertheless, 
during the estimated times of major cladogenesis events 
within the genus (ca. 8—3.8 million years ago (MYA); Ma- 
kokha et al. 2007; Garcia-Porta et al. 2019) these poten- 
tial physical barriers of gene flow may have been less po- 
rous and only later, in concert with the Plio-Pleistocene/ 
Holocene climatic changes, may have permitted passage. 

During the Early Miocene conditions in southern 
Africa were wetter and warmer than today. Since then 
southern Africa has undergone stepwise cooling and 
drying. Diekmann et al. (2003) suggested that arid 
periods occurred in southern Africa between 11.6 and 
10.7 MYA and between 9.7 and 7.7 MYA. The region 
experienced a third major cooling step with increased 
aridification at the end of the Pliocene (around 2.8 
MYA) (Van Zinderen Bakker 1986; Partridge 1998; 
Diekmann 2003). This coincided with a major uplift- 
ing of the subcontinent (the African superswell; Lith- 
gow-Bertelloni and Silver 1998) resulting in major 
changes of erosional patterns due to rejuvenation and 
reactivation of drainage systems and increased fluvial 
activity (Partridge et al. 2006). These drastic climat- 
ic changes may have resulted in newly available open 
habitats for Pedioplanis species occupying more arid 
habitats like P. branchi, facilitating subsequent ra- 
diation and then likely led to range contractions and 
fragmentation among ancestral Pedioplanis spp. with 
subsequent cladogenesis, scenarios that have been sug- 
gested for e.g. southern African Pachydactylus geckos 
(Bauer 1999; Heinicke et al. 2017) and lacertids gener- 
ally (Garcia-Porta et al. 2019). 

The major uplift together with changes in drainage 
systems during less arid periods at the end of the Pliocene 
seems to coincide well with the split of P. inornata and P. 
undata, indicating that the now ephemeral Kuiseb and/or 
Swakop River might have been effective impediments to 
gene flow between these two sister taxa in the past. To- 
day, in the east there are populations of P. undata occur- 
ring south of the mostly dry Kuiseb, indicating that gene 
flow between the two taxa is possible. For P. mayeri (the 
basal taxon of the clade) and P. gaerdesi, the escarpment 
has been suggested to form a geographical barrier (Bau- 
er et al. 1993). During our study we detected a narrow 
area of potential sympatry north of Palmwag and — based 
on museum specimens — near Sesfontein, suggesting that 
while the escarpment might have formed a barrier in the 
past it no longer does. 
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Evolution of body coloration and patterning 
within the P undata species complex 


The P. undata species complex has historically proven 
taxonomically challenging, in part due to the fact that 
most species have been distinguished chiefly by their 
body coloration and patterning (Branch 1998; Conradie 
et al. 2012). Our phylogenetic results show both that di- 
vergent taxa share similar color and patterning character- 
istics, and that intraspecific variation can be considerable. 

Dorsal striping may be a plesiomorphic condition 
within the P. undata complex given that P. mayeri, the 
sister taxon to the remainder of the complex, as well as 
the Angolan outgroup species P. huntleyi and P. haackei, 
are all typically striped. The absence of dorsal striping in 
P. rubens, P. inornata and P. gaerdesi would suggest that 
stripes were lost in lineages that subsequently invaded 
more open habitats (further discussed below) whereas the 
ancestral, striped condition was retained within P. undata. 

Members of the P. undata species complex are chiefly 
allopatrically distributed in regions of Namibia that are 
ecologically distinct; we suggest that apparent color vari- 
ation both among and within species may be a result of 
differences in habitat and associated selective pressures. 
Uniformly and plainly colored taxa including P. branchi, 
P. gaerdesi, P. inornata and P. rubens all occur in arid re- 
gions where they are camouflaged on bare, sandy or rocky 
substrates that lack dense vegetation (Arino et al. 2012; 
Stone and Thomas 2013). In the case of P. rubens, which 
has a very restricted range in the Waterberg Plateau, in- 
dividuals are reddish and are found exclusively on red 
sandstone, whereas P. gaerdesi and P. branchi are found 
in similar habitats in the pro-Namib and adjacent rocky 
outcrops of the Namib, respectively, and each possess- 
es a light gray dorsum, especially anteriorly. Likewise, 
striped patterning in P. mayeri and northern populations 
of P. undata may also be a consequence of shared similar- 
ities in habitat type. Both taxa occur in inland localities in 
northern Namibia, a region that recetves 200-500 mm of 
rainfall annually (Rohde and Hoffman 2012) and in hab- 
itats predominantly characterized by annual and perenni- 
al grasses (Arino et al. 2012; Wang et al. 2012). Among 
squamates the correlation between dorsal striping and the 
use of grassy substrates has been hypothesized to reflect 
the use of disruptive coloration in order to reduce the 
likelihood of detection by predators during escape (Wolf 
and Werner 1994; Farallo and Forstner 2012). 

This hypothesized link between ecology and body col- 
oration in the P. undata species complex may be further 
demonstrated by the pronounced intraspecific variation ev- 
ident within the sister taxa P. inornata and P. undata. For 
example, within P. inornata, typical specimens possess a 
predominantly gray dorsum, however we encountered lo- 
calized red morphs on red, hard-packed soils ca. 60 km 
west of Mariental, and dark morphs on basalt rock piles 
just north of the same town. Pedioplanis undata collected 
in northern localities that are characterized by substantial 
grass cover (Harris et al. 2014) exhibit bold dorsal stripes 
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that appear nearly indistinguishable from those exhibited 
by P. mayeri (see Fig. 4 for images of dorsal color pattern- 
ing typical of the species), whereas individuals collected 
Just south of the Khomas Hochland, where grasses are more 
sparsely distributed, possess faint gray dorsal stripes; fur- 
thermore, individuals collected near the Gobabeb Training 
and Research Centre on vegetationless, schist rock outcrops 
had almost patternless, nearly uniformly gray dorsums. 

Taken together, our phylogenetic results and field ob- 
servations strongly suggest that ecology is a significant 
evolutionary driver of body coloration and patterning 
within the P. undata species complex, as has also been 
demonstrated previously for e.g. P. /ineoocellata (Branch 
1998; Edwards 2013; Childers 2015). 


Conclusion 


The results of our study indicate that species diversity 
within Namibian Pedioplanis is greater than previous- 
ly thought and that this diversity is concentrated within 
the P. undata species complex. These results settle the 
long-standing hypothesis originally put forth by Mayer 
and Berger-Dell’mour (1987) on the existence of addition- 
al lineages within the complex and confirm that 1t compris- 
eS SIx species rather than the four traditionally considered. 

Based on present species distributions and divergence 
date estimates derived from previous phylogenetic studies 
on African lacertids (Makokha et al. 2007; Garcia-Porta et 
al. 2019 we propose several biogeographic scenarios that 
may have led to cladogenesis in the P. undata species com- 
plex. We hypothesize that species diversification beginning 
at ca. 8—3.8 MYA may have been precipitated by major 
climatic changes in combination with the formation of bio- 
geographic barriers across the subcontinent, causing once 
widespread populations to become isolated and restricting 
gene flow until cladogenesis ultimately ensued. Among 
these isolated populations, differences in habitat and as- 
sociated natural selection pressures may have produced 
phenotypic differences among them including the evident 
diversity in body color and patterning between species. 

This study highlights the need for dense geographic 
sampling among broadly distributed groups that have ra- 
diated recently and rapidly, as is the case for the genus 
Pedioplanis. Furthermore, without fine-scaled microhab- 
itat data it is difficult to accurately assess ecological dif- 
ferences and to test our hypothesis of convergent evolu- 
tion in this group, and we recommend future researchers 
prioritize these critical data. 

The majority of organismal research in Namibia has 
so far mainly focused on the Namib Desert region in the 
west (Herrmann and Branch 2013), where species discov- 
ery has accelerated over the past decade due in part to the 
implementation of molecular phylogenetics. However, 
this study highlights the need for increased biodiversity 
surveys of the central interior portion of Namibia, as our 
results suggest that squamate endemicity in this region 
may have been previously underestimated. 
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M350421, M350427, M350432, M350437, M350442). 
Type specimens are denoted with asterisks (** = holo- 
type; * = paratype). All raw data is available upon request 
from the corresponding author (JLC). 

Pedioplanis branchi sp. nov. — CAS 1262218, 1734233, 
1734428, 1734458, 173446, 2147875, 214788'**, 214789, 
2147908* , 2147928* 2147938 2147948*; MCZ-R1841625 
(x-ray only); ZMB 89302, 89303, 89304, 893058*, 89306, 
89307, 89308, 89309, 893108*, 893118*, 89312, 89313, 
89314, 89315, 893168*, 89317, 89318, 89319, 89320, 
89321, 89613. 

Pedioplanis gaerdesi — CAS 206964, 214599, 214729, 
214730, 214741, 214742, 214743, 214744, 214745, 
214746; MCZ-R184294, 184959, 185875, 185883, 
185885, 185887, 185888, 188237; ZMB 89328, 89329, 
89330, 89331, 89332, 89612. 

Pedioplanis inornata —CAS 193404, 200012, 200039, 
200096, 201873, 203522; MCZ-R183772, 184767, 
185890, 193103, 193104, 193105, 193293, 193311, 
193370; ZMB 89333, 89334, 89336, 89337, 89338, 
89339, 89340, 89341. 

Pedioplanis mayeri sp. nov.— CAS 214643'*, 2146448, 
2146458* , 2240108, 224013; MCZ-R1858705*, 185872'*, 
1882238, 1902135*, 1931258*, 1931788, 1931798**; ZMB 
80391, 893498* | 893508*, 89624, 89625. 

Pedioplanis rubens — ZMB 84748, 89351, 89615, 
89616, 89617, 89618, 89619, 89620, 89621, 89622, 
89623. 

Pedioplanis undata — MCZ-R188247, 188248, 
193272, 193276, 193376; ZMB 91631, 89626, 89627, 
89628, 89629, 89630, 89631. 
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Sampling list for all specimens used in 
phylogenetic analyses 


Authors: Jackie L. Childers, Sebastian Kirchhof, Aaron 
M. Bauer 

Data type: excel table 

Explanation note: Locality information and GenBank ac- 
cession numbers are provided for all tissue vouchers 
used in this study. We also provide the specimen ID 
numbers and mtDNA clade assignments for members of 
the Pedioplanis undata species complex that are depict- 
ed in Fig. 2a, b. Institutional abbreviations for catalog 
numbers are as follows: CAS = California Academy 
of Sciences, MCZ = Museum of Comparative Zoolo- 
gy, NHMW = Naturhistorisches Museum Wien, PEM 
= Bayworld (Port Elizabeth, SA), ZMB = Museum fir 
Naturkunde. Collector and associated field number ab- 
breviations are as follows: ABA-ABD-ABE = Werner 
Mayer, AMB = Aaron M. Bauer, JM = Johan Marais, 
JVV = Jens V. Vindum, KTH = Krystal Tolley, MB- 
MBUR = Marius Burger, MH = Michael Cunningham, 
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N = Anamarija Zagar, RBH = Raymond B. Huey, SK 
= Sebastian Kirchhof, SVV = Stuart V. Nielson, USH 
= University of Stellenbosch Herpetology, WC = Wer- 
ner Conradie. AMB, MCZ-A, and MCZ-Z specimens 
lacking institutional registration numbers are currently 
awaiting cataloging at the National Museum of Namibia. 
Copyright notice: This dataset 1s made available under 
the Open Database License (http://opendatacommons. 
org/licenses/odbl/1.0/). The Open Database License 
(ODbL) is a license agreement intended to allow us- 
ers to freely share, modify, and use this Dataset while 
maintaining this same freedom for others, provided 
that the original source and author(s) are credited. 
Link: https://doi.org/10.3897/zse.97.61351.suppl1 
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Mensural data for P. undata species complex 
specimens examined. 


Authors: Jackie L. Childers, Sebastian Kirchhof, Aaron 
M. Bauer 

Data type: excel table 

Explanation note: Minimum and maximum values (mm) 
for each species are reported for each character, see 
Materials and methods for character abbreviations and 
descriptions. Only adults were examined, and only 
individuals with intact tails were used to obtain tail 
length values (TaL); see Appendix 1 for a list of all 
specimens examined. 

Copyright notice: This dataset is made available under 
the Open Database License (http://opendatacommons. 
org/licenses/odbl/1.0/). The Open Database License 
(ODbL) is a license agreement intended to allow us- 
ers to freely share, modify, and use this Dataset while 
maintaining this same freedom for others, provided 
that the original source and author(s) are credited. 

Link: https://doi.org/10.3897/zse.97.61351.suppl2 
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Meristic and scalation data for P undata 
species complex specimens examined. 


Authors: Jackie L. Childers, Sebastian Kirchhof, Aaron 
M. Bauer 

Data type: excel table 

Explanation note: Abbreviations for character values are 
as follows: C = In Contact; S = Separated by granule; 
SF = Separated by frontal and frontonasal; SP = Sep- 
arated by parietals. Data were collected from juvenile 
and adult specimens. See Materials and methods for 
character abbreviations and descriptions; see Appen- 
dix 1 for a list of all specimens examined. 
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Copyright notice: This dataset is made available under 
the Open Database License (http://opendatacommons. 
org/licenses/odbl/1.0/). The Open Database License 
(ODbL) is a license agreement intended to allow us- 
ers to freely share, modify, and use this Dataset while 
maintaining this same freedom for others, provided 
that the original source and author(s) are credited. 

Link: https://do1.org/10.3897/zse.97.61351.suppl3 
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Table S1. Primer pairs 


Authors: Jackie L. Childers, Sebastian Kirchhof, Aaron 
M. Bauer 

Data type: excel table 

Explanation note: Primers pairs used in this study; PCR 
annealing temperatures (°C) are provided. 

Copyright notice: This dataset is made available under 
the Open Database License (http://opendatacommons. 
org/licenses/odbl/1.0/). The Open Database License 
(ODbL) 1s a license agreement intended to allow us- 
ers to freely share, modify, and use this Dataset while 
maintaining this same freedom for others, provided 
that the original source and author(s) are credited. 

Link: https://do1.org/10.3897/zse.97.61351.suppl4 
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Table S2. Bayesian partitioning schemes 


Author: Jackie L. Childers 

Data type: excel table 

Explanation note: Bayesian partitioning schemes for the 
concatenated mtDNA (Partition nos. 1—4) and concat- 
enated nDNA (Partition nos. 5—8) datasets based on 
the results of our PartitionFinder v.2.1.1 analyses, in 
which best-fitttng models of evolution were selected 
using the Bayesian Information Criterion (BIC) score. 

Copyright notice: This dataset is made available under 
the Open Database License (http://opendatacommons. 
org/licenses/odbl/1.0/). The Open Database License 
(ODbL) is a license agreement intended to allow us- 
ers to freely share, modify, and use this Dataset while 
maintaining this same freedom for others, provided 
that the original source and author(s) are credited. 

Link: https://do1.org/10.3897/zse.97.61351.suppI5 
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Table S3. Summary statistics for select 
meristic variables for PR undata species 
complex 


Author: Sebastian Kirchhof 

Data type: excel table 

Explanation note: Summarized results of Shapiro-Wilk 
normality test, Levene’s test, Kruskal-Wallis rank sum 
test and ANOVA for selected meristic variables for 
adults (SVL > 40 mm) of all species of the P. undata 
species complex. Columns show character (see Mate- 
rials and Methods for abbreviations and descriptions), 
Shapiro-Wilk W (Shapiro _W) and _ corresponding 
p-value (Shapiro_p), Levene’s F (Levene_F) and cor- 
responding p-value (Levene_p), Kruskal-Wallis chi- 
square statistics (KW_ chi-squared) and corresponding 
p-value (KW_p), ANOVA F (ANOVA F) and corre- 
sponding p-value (ANOVA_p) and degrees of freedom 
(df). Differences are considered significant with p < 
0.05 and significance codes are as follows: * =p < 0.5; 
B= Oils te" 1p =.0:0011 

Copyright notice: This dataset is made available under 
the Open Database License (http://opendatacommons. 
org/licenses/odbl/1.0/). The Open Database License 
(ODbL) is a license agreement intended to allow us- 
ers to freely share, modify, and use this Dataset while 
maintaining this same freedom for others, provided 
that the original source and author(s) are credited. 

Link: https://do1.org/10.3897/zse.97.61351.suppl6 
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Table S4. Summarized results of Wilcoxon 
rank sum tests for the meristic characters 


Author: Sebastian Kirchhof 

Data type: excel table 

Explanation note: Summarized results of Wilcoxon rank 
sum tests for the meristic characters with significant 
Kruskal-Wallis test results (number of subdigital 
lamellae SuL and number of femoral pores on right leg 
Fe). All possible species pairs (species pair) for spe- 
cies of the P. undata species complex were tested. Col- 
umns show Wilcoxon W (Wilcox_W) and correspond- 
ing p-value (Wilcox_p) and significant differences (p 
< 0.05) are indicated with * and significance codes are 
as follows: * =p <0.5; ** =p <0.01; *** =p<0.001. 

Copyright notice: This dataset is made available under 
the Open Database License (http://opendatacommons. 
org/licenses/odbl/1.0/). The Open Database License 
(ODbL) is a license agreement intended to allow us- 
ers to freely share, modify, and use this Dataset while 
maintaining this same freedom for others, provided 
that the original source and author(s) are credited. 

Link: https://do1.org/10.3897/zse.97.61351.suppl7 
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Figure S1. Linear relationships with 95% 
confidence level intervals of mensural 
characters between P. mayeri (= P. undata 
"North") and P. undata 


Author: Sebastian Kirchhof 

Data type: PNG image 

Explanation note: Linear relationships with 95% con- 
fidence level intervals of mensural characters with 
snout—vent length (SVL) for adult (SVL > 40 mm) P. 
mayeri (= P. undata “North” sensu Mayer and Berg- 
er-Dell’mour 1987) and P. undata. Linear regression 
equations are provided when the slope coefficient was 
significant (p < 0.05), in case of non-significance this 
is indicated by “not sig.” abbreviations; descriptions 
of characters are provided in Materials and methods. 

Copyright notice: This dataset is made available under 
the Open Database License (http://opendatacommons. 
org/licenses/odbl/1.0/). The Open Database License 
(ODbL) is a license agreement intended to allow us- 
ers to freely share, modify, and use this Dataset while 
maintaining this same freedom for others, provided 
that the original source and author(s) are credited. 

Link: https://do1.org/10.3897/zse.97.61351.suppl8 
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Figure S2. Linear relationships with 95% 
confidence level intervals of mensural 
characters between P. branchi (= P. inornata 
Central") and P inornata 


Author: Sebastian Kirchhof 

Data type: PNG image 

Explanation note: Linear relationships with 95% con- 
fidence level intervals of mensural characters with 
snout—vent length (SVL) for adult (SVL = 40 mm) P. 
branchi (= P. inornata “Central” sensu Makokha et al. 
2007) and P. inornata. Linear regression equations are 
provided when the slope coefficient was significant (p 
< 0.05), in case of non-significance this is indicated 
with “not sig.” abbreviations; descriptions of charac- 
ters are provided in Materials and methods. 

Copyright notice: This dataset is made available under 
the Open Database License (http://opendatacommons. 
org/licenses/odbl/1.0/). The Open Database License 
(ODbL) is a license agreement intended to allow us- 
ers to freely share, modify, and use this Dataset while 
maintaining this same freedom for others, provided 
that the original source and author(s) are credited. 

Link: https://doi.org/10.3897/zse.97.61351.suppl9 


